Differential isospin-fractionation in dilute asymmetric nuclear matter 
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The differential isospin-fractionation (IsoF) during the liquid-gas phase transition in dilute asym- 
metric nuclear matter is studied as a function of nucleon momentum. Within a self-consistent 
thermal model it is shown that the neutron/proton ratio of the gas phase becomes smaller than 
that of the liquid phase for energetic nucleons, although the gas phase is overall more neutron-rich. 
Clear indications of the differential IsoF consistent with the thermal model predictions are demon- 
strated within a transport model for heavy-ion reactions. Future comparisons with experimental 
data will allow us to extract critical information about the momentum dependence of the isovector 
strong interaction. 

PACS numbers: 21.65.+f, 21.30.Fe, 24.10.Pa, 64.10.+h 



Properties of isospin asymmetric nuclear matter play 
an important role in understanding many key issues in 
both nuclear physics and astrophysics [l], [3, [1,111 • To ex- 
plore these properties has been one of the major objec- 
tives of nuclear sciences. Recent progress and new plans 
in conducting experiments with more advanced radioac- 
tive beams provided us a great opportunity to achieve 
this goal. Meanwhile, impressive advances have also been 
made recently in several theoretical fronts in exploring 
the nature of neutron-rich nucleonic matter and their ex- 
perimental manifestations. One of the especially interest- 
ing new features of a dilute asymmetric nuclear matter 
is the isospin-fractionation (IsoF) during the liquid-gas 
(LG) phase transition in it @, 0, [Hi HH . The non-equal 
partition of the system's isospin asymmetry with the gas 
phase being more neutron-rich than the liquid phase has 
been found as a general phenomenon using essentially all 
thermodynamical models and in simulations of heavy-ion 
reactions, for reviews see, e.g., Refs. [H IS IS Ell Ell Ell • 

In fact, indications of the IsoF have been reported ever 
since the early 80 's although their interpretations have 
not always been unique [ill, EH- Several recent exper- 
iments confirmed unambiguously the IsoF phenomenon 
@- lu particular, the experiments and analyses by Xu 



et al. l_5j at the NSCL/MSU are among the most inter- 
esting and detailed ones available so far. In their experi- 
ments the isotope, isotone and isobar ratios were utilized 
to obtain an estimate of the neutron/proton density ra- 
tio Pn/Pp in the gas phase at the breakup stage of the 
reaction. It was found clearly that the gas phase was 
significantly enriched in neutrons relative to the liquid 
phase represented by bound nuclei. However, in all of 
the existing theoretical and experimental studies in the 
literature, only the average neutron/proton ratios inte- 
grated over the nucleon momentum in the liquid and 
gas phases were studied. We normally refer the above 
predicted and observed isospin-fractionation as the inte- 
grated IsoF. In this work, we investigate the differential 
IsoF as a function of nucleon momentum. Surprisingly, 
completely new and very interesting physics is revealed 



from the fine structure of the differential IsoF analysis. 
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FIG. 1: (Color online) The symmetry potential and energy 
(insert) in the MDI interaction with x = and x — — 1. 

Our study is carried out first within a self-consistent 
thermal model [l6| using the isospin and momentum- 
dependent MDI interaction [T3, fig . With this inter- 
action, the potential energy density V(p, T, S) at to- 
tal density p, temperature T and isospin asymmetry 
5 = (p n - pp)/p is 

V[p,T,S) = — + — ( Pn +p p ) + _ | - a (1) 
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In the mean field approximation, Eq. |T]) leads to the 
following single particle potential for a nucleon with mo- 
mentum p and isospin r 

U T (p, T, S, p, x) = A u (x) + Ai(x) — 
Po Po 

I (1 - x£ 2 ) - 8TZ 5p- T 

\PoJ cr + 1 pi 



t — T,—T 



Po 
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where r = 1/2 (—1/2) for neutrons (protons), x, A u (x), 
Ai(x), B, C T .. T ,C T ._ T , <r, and A are all parameters given 
in Ref. [171 ] . This interaction gives an incompress- 
ibility of K$ = 211 MeV for symmetric nuclear mat- 
ter at saturation density po = 0.16 fm~ 3 . The dif- 
ferent x values in the MDI interaction are introduced 
to vary the density dependence of the nuclear symme- 
try energy while keeping oth er p roperties of the nuclear 
equation of state fixed [l7], [HI]- It is worth mention- 
ing that the nucleon isoscalar potential estimated from 
EAsoscaiar ~ (U n + U p )/2 agrees with the prediction of 
variational many-body calculations for symmetric nu- 
clear matter [19| in a broad density and momentum range 
[l7l ] . The corresponding isovector (symmetry) potential 



can be estimated from U s 



(U n -U p )/25. With both 



x = and x = — 1 at normal nuclear matter density, the 
symmetry potential agrees very well with the Lane po- 
tential extracted from nucleon-nucleus and (n,p) charge 
exchange reactions available for nucleon kinetic energies 
up to about 100 MeV [13, H3]- At abnormal densities and 
higher nucleon energies, however, there is no experimen- 
tal constrain available at present. 

The isospin diffusion data from the NSCL/MSU 
have recently constrained the value of x to be between 
and —1 at sub-saturation densities [IHI12]- Shown in 
the insert of Fig. Q] is the density dependent symme- 
try energy E sym (p) with x — and x = — 1, respectively. 
While this constraint on the E sym (p) is the most stringent 
achieved so far in the field, the corresponding symmetry 
potential shown in Fig. [T] still diverges much more widely 
with both momentum and density. This is not surprising. 
Comparing Eqs. |T]) and ([2]), it is seen that the symmetry 
energy obtained using Eq. ([1]) involves the integration of 
the single nucleon potential over its momentum. To ob- 
tain information about the underlying momentum- and 
density-dependence of the symmetry potential which is 
more fundamental than the E sym (p) for many important 
physics questions, one thus has to use differential probes. 
We will demonstrate that the differential IsoF as a func- 
tion of nucleon momentum is such an observable. 

The phase space distribution function f T is the Fermi 
distribution 

2 

fr(r,p) = ^[exp(^^ ) + l]- 1 (3) 

where p T is the proton or neutron chemical potential de- 
termined self-consistently from p T = J f T (r,p)d 3 p. In the 
above, m T is the proton or neutron mass. From a self- 
consistency iteration scheme [HI US], the chemical po- 
tential \i T and the distribution function f T (f,p) are de- 
termined numerically. 

From the chemical potential p T and the distribution 
function f T (r,p), the energy per nucleon E(p,T,S) can 
be obtained as 



E(p,T,5) = - 
P 



Furthermore, the entropy per nucleon S T {p,T 1 8) is ob- 
tained from 



S T {p,T, 5) 
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Jo 

(5) 

with the occupation probability n T = 4-/t(^*)P)- Finally, 
the pressure P(p, T, S) is calculated from the thermody- 
namic relation 



P(p,T,S) 



Tj2S T (p,T,5)-E(p,T,5) 



PtPt 



(6) 

For the ease of following discussions, we first reproduce 
the integrated IsoF using the above formalism. In hot 
asymmetric nuclear matter, the LG phase coexistence is 
governed by the Gibbs conditions 



^(T,pf) 
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At any given pressure below the critical point, the two so- 
lutions for the liquid and gas phases thus form the edges 
of a rectangle in the proton and neutron chemical poten- 
tial isobars as a function of isospin asymmetry 6 and can 
be found by means of the standard geometrical construc- 
tion method fijl . 
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FIG. 2: (Color online) Left window: the section of binodal 
surface at T = 10 MeV with x = and x = -1. The crit- 
ical point (CP), the points of equal concentration (EC) and 
maximal asymmetry (MA) are also indicated. Right window: 
the double neutron/proton ratio in the gas and liquid phases 
{n/p)a I '(n/p)L as a function of the nucleon kinetic energy. 

Shown in the left window of Fig. [2] is a typical section 
of the binodal surface at T = 10 MeV with x = and 
x = — 1. The phenomenon of integrated IsoF with the 
gas phase being more neutron-rich is seen very clearly. It 
is also seen that the stiffer symmetry energy {x = — 1) 
significantly lowers the critical point (CP). However, be- 
low a pressure of about P = 0.12 MeV/fm 3 , the magni- 
tude of the integrated IsoF becomes almost independent 
of the symmetry energy used. To examine the advan- 
tages of the differential IsoF analyses over the integrated 
ones, as an example, we select the gas and liquid phases 



3 



in equilibrium at T=10 MeV and P = 0.1 MeV/fm 3 . 
For x — the density and isospin asymmetry are, re- 
spectively, pa — 0.087po and 8 G — 0.791 for the gas 
phase, and pz = 0.763/?o and 5l — 0.296 for the liquid 
phase. For x = — 1 they are, respectively, p G = 0. 114/30, 
S G = 0.808, pl = 0.714p and S L = 0.30. The cor- 
responding double neutron/proton ratio in the gas and 

liquid phases (^/p)^ (p) can then be readily obtained us- 
ing Eq. ([3]) as a function of nucleon momentum or kinetic 
energy. We refer this function as the differential IsoF. 

Shown in the right window of Fig. [5] are the differen- 
tial IsoFs for both x = and x = — 1. It is interesting to 
see that the isospin-fractionation is strongly momentum 
dependent. Moreover, while the integrated double neu- 
tron/proton ratios of 5.07 (x = —1) and 4.65 (x = 0) are 
very close to each other, the differential IsoF for nucle- 
ons with kinetic energies high than about 50 MeV is very 
sensitive to the parameter x used. Surprisingly, a rever- 
sal of the normal IsoF is seen for x = for nucleons with 
kinetic energies higher than about 220 MeV. In this case, 
there are more energetic neutrons than protons in the liq- 
uid phase compared to the gas phase. We note that at 
pressures higher than 0.1 MeV/fm 3 where the integrated 
IsoF is already very sensitive to the E sym (p) , the differen- 
tial IsoF is much more sensitive to the x parameter than 
that shown in Fig. O For energetic nucleons where the 
differential IsoF is very sensitive to the parameter x, the 
f T can be well approximated by the Boltzmann distribu- 
tion as shown in Fig. [5] For these nucleons in either the 
liquid (L) or gas (G) phase, the neutron/proton ratio 

(n/p) L/G = eM^ /G ~E^ G -p^ G +p L p ' G )/T]. (9) 

The energy difference of neutrons and protons having the 
same kinetic energy and mass 

E L/G _ E L/G = jjL/G _ jjL/G K 2( j . Usym (j>, p L/G ) 

(10) 

is directly related to the symmetry potential U sym . Be- 
cause of the chemical equilibrium conditions given in Eq. 
([7]) , the chemical potentials cancel out in the double neu- 
tron/proton ratio 

G {p) = c-xp[-2(d G -U S y m (p,pG)-SL-Usy m (p, Pl))/T] 

(11) 

This general expression clearly demonstrates that the dif- 
ferential IsoF for energetic nucleons carries direct infor- 
mation about the momentum dependence of the symme- 
try potential. In the above expressions, the tiny tem- 
perature dependence of the symmetry potential has been 
neglected. 

For both the liquid-gas and hadron-QGP (quark- 
gluon-plasma) phase transitions, equilibrium model cal- 
culations for infinite nuclear matter are very useful for de- 
veloping new concepts and predicting novel phenomena. 
However, the experimental search/confirmation for the 
new phenomea/concepts in real nuclear reactions is usu- 
ally very challenging. For example, the intermediate en- 
ergy heavy-ion reaction community has been studying for 



{n/p)i 



more than two decades the underlying nature and exper- 
imental signatures of the LG phase transition predicted 
first for infinite nuclear matter based on thermodynami- 
cal considerations. Hopefully, similar to the very fruitful 
study recently in momentum-space about how electrons 
behave during t he p hase transition to high-temperature 
superconductors [24j, the new study about how nucleons 
behave in the correlated momentum- and isospin-space 
may reveal deeper insights into the nature of the LG 
phase transition. 
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FIG. 3: (Color online) The neutron/proton ratio in the "gas" 
and "liquid" phases as a function of the nucleon kinetic energy 
for the reaction of 124 Sn+ 124 Sn at E beam /A = 50 MeV. 
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FIG. 4: (Color online) Same as Fig. 3 but for the reaction of 
i32 Sn+ i24 Sn at Ebeam / A = 400 MeV. 



How can one experimentally measure the differential 
IsoF? As for the structure functions of quarks and glu- 
ons in the initial state of relativistic heavy-ion collisions, 
the momentum distribution of the n/p ratio in the liquid 
phase may not be measured directly since what can be de- 
tected at the end of heavy-ion reactions are free nucleons 
and bound nuclei in their ground states. Nevertheless, 
precursors and/or residues of the transition in the differ- 
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ential IsoF may still be detectable in heavy-ion reactions 
especially those induced by radioactive beams. While a 
comprehensive study on the detailed signatures of the dif- 
ferential IsoF is beyond the scope of this work, we have 
actually performed calculations using the IBUU04 trans- 
port model with the same MDI interaction!!?]]- Shown 
in Fig. [3] and Fig. 2] are two typical examples for the 
central reaction of 124 Sn+ 124 Sn at E beam /A = 50 MeV 
and 132 Sn+ 124 Sn at E beam /A = 400 MeV, respectively 
To separate approximately nuclcons in the low density 
"gas" region from those in the "liquid" region a density 
cut at 0.13yOo is used. Calculations using a cut at 0.5po 
in the reaction at 50 MeV/ A lead to qualitatively the 
same although quantitatively slightly different results. 
Very interestingly, in both reactions there is indeed a 
transition from the neutron-richer (poorer) "gas (liquid)" 
phase normally known as the IsoF for low energy nucleons 
to the opposite behavior (i.e., anti-IsoF) for more ener- 
getic ones. Moreover, the transition nucleon energy from 
the normal IsoF to the anti-IsoF is sensitive to the pa- 
rameter x used. This is more pronounced in the reaction 
at E}, eam /A — 400 MeV where effects of the symmetry 
(Coulomb) potential are relatively stronger (weaker) for 
more energetic nucleons consistent with predictions of the 
thermal model in the right window of Fig. 2. Comparing 
the thermal model predictions and the transport model 
results, one sees that the two independent approaches 
predict qualitatively the same phenomenon while there 
are quantitative differences, especially for low energy nu- 
cleons. This is mainly because in nuclear reactions the 
Coulomb repulsion shifts protons in the gas phase from 
low to higher energies leading to the peak in (n/p)c ratio 
at Ekin — 0, while it has little effects on the protons in 
the liquid phase. Transport model calculations without 
the Coulomb potential lead to results more closely resem- 
bling the thermal model predictions for infinite nuclear 
matter. We notice that the "gas" phase defined in this 



study contains also the pre-equilibrium nucleons which 
are known to be more neutron-rich than the reaction 
system 125]. They are energetic and thus affect mostly 
the high energy part of the (n/p)Q ratio. The subtrac- 
tion of the pre-equilibrium nucleons from our analyses 
will thus mainly lower the {ji/p)g for high energy nuclc- 
ons. Therefore, the transition from the normal IsoF to 
the anti-IsoF will become more obvious and it will make 
our conclusions even stronger. We also notice that the 
radial flow affects all nucleons in both the gas and liquid 
phases the same way as already illustrated in ref . [2q . 
The subtraction of the radial flow is expected to have 
very little effect on the neutron/proton ratio. 

In summary, the differential IsoF is studied as a func- 
tion of nucleon momentum. While the gas phase is overall 
more neutron-rich than the liquid phase, the gas phase is 
richer (poorer) only in low (high) energy neutrons than 
the liquid phase. Clear indications of the differential IsoF 
consistent with the thermodynamic model predictions are 
demonstrated within a transport model for heavy-ion re- 
actions. While it may be very challenging to test experi- 
mentally the predictions, future comparisons with exper- 
imental data will allow us to extract critical information 
about the momentum dependence of the isovector strong 
interaction. 
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